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Abstract 

The problem we consider in this work is to minimize the L'^-norm 
{q > 2) of the semidiscrete controls. As shown in |LT06) . under the 
main approximation assumptions that the discretized semigroup is 
uniformly analytic and that the degree of unboundedness of control 
operator is lower than 1/2, the uniform controllability property of 
semidiscrete approximations for the parabolic systems is achieved in 
L^. In the present paper, we show that the uniform controllability 
property still continue to be asserted in L'^{q > 2) even with the con- 
dition that the degree of unboundedness of control operator is greater 
than 1/2. Moreover, the minimization procedure to compute the ap- 
proximation controls is provided. An example of application is imple- 
mented for the one dimensional heat equation with Dirichlet boundary 
control. 

1 Introduction 

Consider an infinite dimensional linear control system 

yit) = Ayit) + Buit), y(0) = yo, (1) 

where the state y(t) belongs to a reflexive Banach space X, the control u(t) 
belongs to a reflexive Banach space U, A : D{A) — )■ X is an operator, and 
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B is a control operator (in general, unbounded) on U. Discretizing this par- 
tial differential equation by using, for instance, a finite difference or a finite 
element scheme, leads to a family of finite dimensional linear control systems 



yh{t) = AnVhit) + BhUhit) 



yh{0) = yoh, 



(2) 



where yhit) G and Uh{t) & Uh, for < h < ho . 

Let yi e X. If the control system (1) is controllable in time T then there 
exists a solution y{.) of (1) associated with a control u such that y(T) = yi. 
As known, we have many available methods in order to identify the controlla- 
bility. We refer to J.-L. Lions |L88j for a well-known method in attainning the 
control of the minimal norm- the so called HUM (HUM stands for Hilbert 
uniqueness method). In this work, however, we investigate a method which 
we can achieve the minimization procedure in L'^ norm (g > 1). Namely, we 
will establish some conditions ensuring the existence and convergence of the 
discretized control of the minimal L'^ norm 



Necessary conditions for optimal control in finite dimensional state space 
were derived by Pontryagin et al |PS62j (see also [E], [T05J). The Maximum 
Principle as a set of necessary conditions for optimal control in infinite di- 
mensional space has been studied by many authors. Since it is well known 
that the Maximum Principle may be false in infinite dimensional space, there 
are still many papers that give some conditions to ensure that the Maximum 
Principle remains true. It was Li and Yao |LY85j who used the Eidelheit sep- 
aration theorem and of the Uhl's theorem in order to extend the Maximum 
Principle to a large class of problems in infinite dimensional spaces when the 
target set is convex and the final time T is fixed. Additionally, the authors 
of [F87j . |FF91j . |LY91] . by making use of Ekeland's variational principle, 
give some conditions on the reachable set and on the target set in order to 
get an extension of Maximum Principle. Nevertheless, the problem is that 
when we applied the result of |F87] . |LY91j for the system (1) in the case 
the final state and final time are fixed, the finite-codimensional condition in 
|F87j . |LY91j does not satisfy for the system (1) in general. Hence we cannot 
use Maximum Principle in our problem. Fortunately, thank to the Fenchel- 
Rockafellar duality theorem which is used in the same manner in |CGL94] . 
|GLH08] . the constrainted minimization of the function can be replaced by 
the unconstrainted minimization problem of corresponding conjugate func- 
tion. Therefore, we will consider the above minimization procedure in the 
same framework with |GLH08] . 




(3) 
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The uniform controllability is an important area of control theory re- 
search and it has been the subject of many papers in recent years. The main 
goal of this article is to establish conditions ensuring a uniform controlla- 
bility property of the family of discretized control systems (2) in L'^ and to 
establish a computationally feasible approximation method for identifying 
controllability. 

It is well known that controUabilty and observability are dual aspects 
of the same problem. We therefore will focus on the uniform observability 
which is shown to hold when the observability constant of the finite dimen- 
sional approximation systems does not depend on h. Some relevant references 
concerning this property has been investigated by many authors in series 
of articles [IZ99] . [LZ98], jLZ02] . [NZ03] . |Zua99j . |Zua02j . |Zua04j . [ZuaOSj . 
|Zua06j . |BHL10a] . |BHL10b] and [LT06] . For finite difference schemes, a uni- 
form observability property holds for one-dimensional heat equation jLZ02] . 
beam equation |LZ98j . Schrodinger equations |Zua05] . but does not hold for 
1-D wave equations |IZ99j . This is due to the fact that the discrete dynam- 
ics generates high frequency spurious solutions for which the group velocity 
vanishes that do not exist at the continuous level. To overcome these high 
frequency spurious for wave equations, [ZuaOS] showed some remedies such 
that Tychonoff s regularization, multigrid method, mixed finite element and 
filtering of high frequency, etc. 

To our knowledge, in 1-D heat equation case, due to the fact that the dis- 
sipative effect of the 1-D heat equation acts as a filtering mechanicsm by itself 
and it is strong enough to exclude high frequency spurious oscillations |LZ98] . 
However, the situation is more complex in multi-dimensional. The counter- 
example is shown in |Zua06] for the simplest finite difference semi-discretization 
scheme for the heat equation in the square. 

In recent works in L^-norm, by means of discrete Carleman inequali- 
ties, the authors in |BHLinaj . |BHLinb] obt ain the weak uniform observ- 
ability inequality for parabolic case by adding reminder terms of the form 
g-c/i ||^/'h(7')||^2(f2) which vanishes asymptotically as — j- 0. Moreover, as 
in |LT06] . the approximate controllability is derived from using semigroup 
arguments and introducing a vanishing term of the form \\4'h{T)\^j^2(p^-^ for 
some /9 > 0. 

In fact, an efficient computing the null control for a numerical approx- 
imation scheme of the heat equation is itself a difficult problem. Accoding 
to the pioneering work of Carthel, Glowinski and Lions in |CGL94] . the null 
control problem is reduced to the minimization of a dual conjugate function 
with respect to final condition of the adjoint state. However, as a consequence 
of high regularizing property of the heat kenel, this final condition does not 
belong to L^, but a much large space that can hardly be approximated by 
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standard techniques in numerical analysic. Recently, A. Munch and collabo- 
rators have developed some feasible numericals such that the transmutation 
method, variational approach, dual and primal algorithms allow to more ef- 
ficiently compute the null control (see in series |CM09] . [CMIO] , |MZ10] . 

IPMinp . 

The discretization framework in this paper is the same spirit as |LT06] . 
|LT00j ■ In |LT06j . under standard assumptions on the discretization process 
and for an exactly null controllable parabolic system (1), if the degree of 
unboundedness of the control operator is lower than 1/2 then the 
semidiscrete approximation models are uniformly controllable and they also 
showed that for the (2), the minimizing of the cost function of discretized 
control with power g = 2 is obtained. 

In this article, we prove the existence of the minimum of the cost function 
of discretized control power q {q > 2) for type (2), in the case the operator A 
generates an analytic semigroup. Of course, due to regularization properties, 
the control system (1) is not exactly controllable in general. Hence, we focus 
on exact null controllability. Our main result, theorem 3.1, states that for 
exactly null controllable parabolic system (1) and under standard approxi- 
mation assumptions, if the discretized semigroup is uniformly analytic, and 
if the degree of unboundedness of the control operator B with re- 
spect to A is greater than 1/2, then a uniform observability inequality in 
(L^) is proved. We stress that we do not prove uniform exact null controlla- 
bility property for the approximating system (2). Moreover, a minimization 
procedure to compute the approximation controls is provided. 

The outline of the paper is as follows. In Section 2, we briefly review some 
well-known facts on controllability of linear partial differential equation in 
reflexive Banach spaces. In Section 3, we consider the existence and unique 
solution of the minimization problem in continuous case. By making use of 
the Fenchel-Rockafellar duality theorem, we gives a constructive way to build 
the control of minimal L'^ norm. The main result is stated in Section 4 and 
proved in Section 5. An example of application and numerical simulations are 
provided in Section 6, for the one-dimensional heat equation with Dirichlet 
boundary control. An Appendix is devoted to the proof of a lemma. 
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2 A short review on controllability of linear 
partial differential equations in reflexive Ba- 
nach spaces 

In this section, it is convenient to first have a quick look to controllabihty 
of infinite dimensional linear control systems in reflexive Banach spaces (see 
more [CTnH] . jP83] . \TWQ7\ ). 

The notation L{E, F) stands for the set of linear continuous mappings 
from E to F, where E and F are reflexive Banach spaces. 

Let X be a reflexive Banach space. Denote by <, >x the inner product on 
X, and by the associated norm. Let S{t) denote a strongly continuous 
semigroup on X, of generator [A, D{A)). Let X_i denote the completion of 
X for norm = — y4)~^a;||, where /3 G p{A) is fixed. Note that X_i 

does not depend on the specific value of /3 G p{A). The space X_i is isomor- 
phic to {D{A*)y, the dual space of D{A*) with respect to the pivot space 
X, and X C X^i, with a continuous and dense embedding. The semigroup 
S{t) extends to a semigroup on still denoted S(t), whose generator is 
an extension of the operator A, still denoted A. With these notations, A is a 
linear operator from X to X_i. 

Let U he a reflexive Banach space. Denote by <,>u the inner product 
on U, and by ||.||^ the associated norm. 

A linear continuous operator B : U ^ X^i is admissible for the semigroup 
S{t) if every solution of 

y' = Ayit) + Bu{t), (4) 

with y{0) = E X and u{.) G /."^(O, +00; f/), satisfies y{t) G X, for every 
t > 0. The solution of equation (1) is understood in the mild sense, i.e, 

y{t) = S{t)y{0) + [ S{t- s)Bu{s)ds, (5) 
Jo 

for every t > 0. 

For T > 0, define Lt : L«(0, T; U) X_i by 

Ltu= [ S{T - s)Bu{s)ds. (6) 
Jo 

A control operator B G L{U, X^i) is admissible, if and only if IttiLt C X, 
for some (and hence for every) T > 0. 
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The adjoint of satisfies 

L*rr:X*-^ (L«(0, T; U))* = 1^(0, T; U*) 
L*T%l,{t) = B*S{T - tyip (7) 
, a.e on [0,T] for every ijj G D{A*). Moreover, we have 

= sup [ < B*S{T - sy^,u{s) >ds, (8) 

for every ip & X* . 

Let B G L{U,X^i) denote an admissible control operator. 

We use two following lemmas (for the proofs we refer to [TW07j ) 

Lemma 1. Z,X are reflexive Banach spaces. G G L{Z,X) then the following 
statements are equivalent: 

• G is onto. 

• G* bounded from below i. e there exists C > such that 

\\G*x\\^ > C \\x\\-^ every x E X. 

Lemma 2. Zi, Z^^ Z^ are reflexive Banach spaces. And f G L{Zi, Z^) , g G 
L{Z2,Z3). Then the following statements are equivalent: 

• Imf C Img. 

• There exists a constant C > such that : \\f*z\\^^ < C\\g*z\\^^^ for 
every z E Z^. 

• There exists an operator h G L{Zi, Z2) such that f = gh. 
We state some concepts as follows 

For j/o G X, and T > 0, the system (4) is exactly controllable from ?/o 
in time T if for every yi G X, there exists u{.) G L'^{0,T]U) so that the 
corresponding solution (4), with y{0) = yo satisfies y(T) = yi. 

In fact that the system (4) is exactly controllable from yo in time T if and 
only if Lx is onto, that is ImLrp = X. Making use of Lemma 1, there exists 
C > such that 

CIlVlIx < II^T^II = / < B*S{T - sy^,u{s) > ds 

< sup [ \\B*S{T -tym\u{t)\\dt 

< { r \\B*s{ty^fdt)p. 

Jo 
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Therefore, the system (4) is exactly controllable in time T if and only if 

r\\B*smrdt>cur^. (9) 

^0 



For T > 0, the system (4) is said to be exactly null controllable in time T 

if for every yo G X, there exists u{.) G L''(0,T; U) so that the corresponding 
solution of ( 4), with y{0) = yo satisfies y{T) = . 

This means that the system (4) is exactly null controllable in time T if and 
only if IniS{T) C ImLx- Making use of Lemma 2 and the same argument 
as above, there exists C > such that 

cwsmw^ < \\L*^^\\ < ( r\\B*s{tr^rdt)-^. 

Jo 

Thus, the system (4) is exactly null controllable in time T if and only if 

r \\B*S{t)*^''dt > C ||5(T)*V^||^ . (10) 
Jo 

3 Duality 

The goal of this section is to show that, with using duality arguments 
and Fenchel- Rockafellar theorem we can achieve the control of minimal 
norm (g > l)for the continuous framework. 

Consider the system : 

'y{t) = Ay{t) + Bu{t) on Qt = (0, T) x Q 

m-yo 

where B is admissible and A generates an analytic semigroup S{t) in the 
reflexive Banach space X. 

Our aim is to mimimize the following functional: 



Minimize J{u) = \ M'' dt (q > 1) 



q J 

Subject to It G -B 

where E = {u eU : u steering the system from y^ at time zero to y(T)=0} 
Theorem 1. The problem (12) has a unique solution u. 
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Proof. First of all, we show the existence of the solution of the optimal prob- 
lem. 

Consider a minimizing sequence {un)neN of controls on [0, T], i.e, 



/ converges to inf J('u) as n ^ +00. (13) 

Jo 

Hence, {un)neN bounded in L^(0,T; U). 

Since U is reflexive space and q < +00, then -L^(0, T; U) is also reflexive. 

Thus, up to a sequence, Un converges weakly to uin L'^. Note that the tra- 
jectory Un (resp. y) associated with the control Un (resp. u) on [0,T] through 
the system 

yn = Ayr, + Bun, i/„(0) = yo, 
and the solution of the above system is expressed in form 

yn{t) = S{t)yo + S{T - s)Bunds. 
Jo 

A passage to the limit imples that 

y = Ay + Bu, 2/(0) = yo, 
and the solution y associated with control u in the form 

y{t)^S{t)yo+ f S{T-s)Buds. 
Jo 

As Un converges weakly to m in L^, we get the inequality 

/ \\u{t)\\Ut < limmi [ \\unit)fdt 
Jo Jq 

= inf I \\u{t)\\'^dt. 
Jo 

It follows easily that 

f \\u{t)\\''dt^\ni r \\u{t)\\Ut. 
Jo Jo 

Hence u is optimal of (12). This ensures the existence of a optimal control. 
Moreover, the cost function is strictly convex then the solution is obvious 
uniqueness. □ 
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By making use of convex duality, the problem of control to trajecto- 
ries is reduced to the minimization of the corresponding conjugate function. 
Roughly speaking, it is stated through the following theorem: 

Theorem 2. (i) We have the identity: 

inf- / \\u\\Ut^-mi{- [ \\B*7j;\fdt+ <jp{0),yo>), (14) 
u&Eq Jq i,T p Jo 

where ip be solution of : 

-ip = A*ip (15) 
m = V'T. (16) 

Or, we have in the form 



inf-/" \\u\\Ut ^ ~ inf (- [ \\B*S(T -tyMPdt+ < S(T)*iP,yo>). 

(17) 

(a) Ifuop is optimal of the problem (12) then 

Uop{t) = \\B*S{T - t) Vf"^ B*S{T - tyip, 
where (p be optimal of the function: 

1 r 

J*{jIj) = - / \\B*S{T-tyiljfdt+ < S{TyiP,yo > . 
P Jo 

Proof. 

(i) Let y be solution of (1) with u = and we introduce the operator 
N e L{Li{Qt),X) with Nu = Zu{.,T) for all u e Li{Qt), where is 
solution to 

z = Az + Bu (18) 
z{x,0) = 0. (19) 

Accordingly, the solution y of (11) can be decomposed in the form 

y^zu + y. (20) 

The adjoint N* is given as follows 

For each i/jt G X*, N*'ipT = B^ip where is solution of (15), (16). 
Let us introduce the following functions F and G 
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1 +00 otherwise ' 

G{u) = - [ \\ufdt. 
q Jo 

Then, the problem (12), where the infimium is taken over all u satisfying 
E, is equivalent to the following minimization problem 

inf (FiNu) + Giu)). (21) 

ueLi{QT) 

We can apply now duality theorem of W.Fenchel and T.R.Rockafellar (see 
Theorem 4.2 p.60 in |ET99] l It gives 

inf {F{Nu) + G{u)) = - inf {G*{N*tPT) + F*{-^t)), (22) 

where F* and G* are the convex conjugate of F and G, respectively. Denote 
that ipT = i^iT), zt = z{T). 
Note that 



F*{^t) = sup < Zt, tpT>= - < ^T, Vt >, 
ZT=-yT 



G*(u) = - r Wcordt. 



for all e X*. 
Additionally, 



P Jo 

Therefore, 

G*(iV>T) + F*(-V^T) = - [ \\B*i;fdt+<M^),yTix)>. (23) 

Finally, multiplying the state equation (15) by y and due to (11), we 
obtain 

< ^T,2/T > = < V'(0),l/0 > • 

Rewrite (23) as follows 

G*{N*tPT) + F*{-ijT) = - r \\B*ijrdt+<ij{0),yo> 

P Jo 

= - [ \\B*SiT-tyM''dt+<SiTyi;T,yo>, 
P Jo 

since ip is the solution of (15) , (16). 
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Prom (21) and (22), we have the identity 



inf- [ Ikll^di = -inf(- /" \\B*M^dt+ <xlj(0),yo>), 
uGEq Jo Vt p Jo 



where ip be solution of (15), (16). 
Or, we have in the form 



inf - [ \\uf dt^- inf (- [ \\B*S(T - tYM^ dt+ < 5(T) Vo >] 



(ii) If we denote by (uop) , ((^t) the unique solution to "LHS of (14)" and 
"RHS of (14)" respectively, then one has 

= -/ \\uop\\'dt+- [ \\B*ifT\fdt+<<fTiO),yo>. (24) 
Q Jo P Jo 

We apply Young's inequality for the first two terms of RHS (24) 

-/ \\uop\\'^dt + - I \\B*^T\fdt> I Uop.B*^T- (25) 

Q Jqt P JQt JQt 

Then, "RHS of (24)" > ^^^Uop.B*VT+ < M^),yo >■ 
Futhermore, by multiplying two sides of (15) by y and applying Green's 
formula, we obtain 

< B*ipT, u> + < (/^r(0), yo >= 0. (26) 

On the one hand, "RHS of (24)" > ( due to (26)). 

On the other hand, "RHS of (24)" = ( due to (24)). 

This is equivalent to that the sign " = " in inequality (14) happens, i.e, 

II Wl II ID* IIP 

\\Uop\\ = \\B (firW . 
It is also rewritten as follows 

Uopit) = \\B*S{T - tYiff'^ B*S{T - tyip, 
where ip be optimal of the function J* is given as above. 

□ 



11 



Remark: It is easily seen that the functional J* is convex, and from the 
inequality (10), is coercive. Then, it follows that J* attains a unique minimum 
in some point (p G D{A*). As above explanation, the control u is chosen by 



u{t) = \\B*S{T - t) Vir B*S{T - t) V, (27) 



for every t G [0,T] and let y(.) be the solution of (11), such that y{0) = yo, 
associated with the control it, then we have y(T)=0. 

Therefore, u is the control of minimal of L'^ norm, among all controls 
whose associated trajectory satisfied y{T) = 0. 

We emphasize that observability in norm (1 < p < 2) implies control- 
lability and gives a constructive way to build the control of minimal L'^ norm 
(g > 2). A similar result was known in norm through using HUM ( see 
more |CT06j . [L88] . |Zua04j . |Zua05j l 

4 The main result 

We are concerned in this work with the uniform controllability property 
for the parabolic systems. As shown in |LT06j . this property is known to 
hold with the degree of unboundedness of control operator 7 G [0, 1/2). In 
this section, we also establish some appropriate assumptions and conditions 
ensuring that the unform controllability still holds in the case 7 G 1/2, 

Let X and U be Hilbert spaces, and let A : D{A) — )■ X be a linear 
operator and self-adjoint, generating a strongly continuous semigroup S{t) 
on X. Let B G L[U, D[A*)') be a control operator. We make the following 
assumptions that will be used along this article (also refer to |LT06j ) 

(HI) The semigroup S(t) is analytic. 

Therefore, (see |P83] ) there exist positive real number Ci and fl such that 

Ult 

\\S{t)\\^ ^ C,e-' \\y\\^ , \\ASit)y\\^ ^ C,— \\y\\^ , (28) 

for alH > and y G D{A), and such that, if we set A = A — ul, for 9 G [0, 1] 
and there holds 

-A')Sit)y ^<C,^\\y\\^, (29) 

for alH > and y G D{A). 

Of course, inequalities (28) hold as well if one replaces Ahj A*, S{t) by 
S{t)*, for y G D{A*). 
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Moreover, if p{A) denotes the resolvent set of A, then there exists S e 
(0, f ) such that p{A)dAs ^ {uj + pe'^\0 > 0, |e| < f + S}. 

For A e p{A), denote by R{\,A) = (A/ - A)~^ the resolvent of A . It 
follows from the previous estimates that exists C2 > such that 



\\R(\mLix) 



< 



c. 



\X — U! 



■,\\AR{X,A)\\ 



L{X)) 



<C2, 



(30) 



for every A e A5, and 

R(\,A)) 



L{X) |A| 



AR{\A) 



L{X) 



<C2, 



(31) 



for every A G + cu. Similarly, inequalities (30) and (31) hold as well with 
A* and A*. 

(H2) The degree of unboundedness of B is 7. Assume that 7 e 1/2^^ 
( where p,q are conjugate, i.e - + - = 1 and 1 < p < 2). This means that 



p q 

B e L{u,D{{-Ayy). 



(32) 



In these conditions, the domain of B* is D{B*) — D{{—A*)'^), and there 
exists C3 > such that 



(33) 



for every ip e D{{-A*)'f). 

(H3) We consider two famihes (-'^/l)o</l</lo {Uh)o<h<ho of finite di- 
mentional spaces, where h is the discretization parameter. 

For every h e (0, ho), there exist the linear mappings Ph : D{{—A*)2y — >■ 

Xh and Ph-. Xh ^ D{{-A*)h) and (i*)"^^^ : D(-(i*)5) ^ D(-(i*)T) 
(resp., there exist linear mappings Qh : U ^ and Qh '■ Uh ^ U), satisfying 
the following requirements: 

(ifs.i) For every h G (0, ho). The following properties hold 

PhPh = idxh and QhQh = iduh- (34) 

(7/3.2) There exist s > and C4 > such that there holds, for every 
h e (0, ho) , 



I -{A* 



(35) 



X 
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for every t/j e D{A*) and 



<W^'-^^\\A*^^, (36) 



^0, 



(J - QhQh)u 
for every u & U, and 

(7 - QhQh)B*^P < C^h<'-^^\\A*i^\\^, 



u 



(37) 



(38) 



for every ip G D(A*) 

For every h G (0, /iq), the vector space (resp. Uh) is endowed with the 
norm ||.||^ (resp. ) defined by: 



PhUh 



Therefore, we have the properties 



for Uh e Xh (resp. ,||Mfe||c;, = QhUh )• 



Qy. 



L{Uh,U) 



1 and 



(i*) ''^^x 



< Clbl 



X 1 



\\Ph\\L[x,Xh) - ^5 and < C5. 

(7^3.3) For every h e (0, h^), there holds 

P/j = Ph and Q;^ = Qh , 



(39) 
(40) 

(41) 



where the adjoint operators are considered with respect to the pivot spaces 
X, U , Xh, Uh- 

(Hsa) There exists Cq such that 



(42) 



for all h G (0, Hq) and iph G Xh- 

For every h G (0, /iq), wc define the approximation operators : Xh 
Xh of A* and B; : Xh ^ Uh of S*, by 



A* = PhA*Ph and = g;,S*(i*) '^''Ph. 



(43) 



(7^4) The following properties hold: 

(7^4.1) The family of operators e*"^^ is uniformly analytic, in sense that 
there exists C7 > such that 



\L{Xh) 
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(44) 



\\^^^'''±ix,)<C7- (45) 

for alH > and h G (0, ho). 

{H42) There exists Cg > such that, for every / G X and every h G 
(0, ho), the respective solutions of A*i/j = f and A'l^ph = Phf satisfy 

\\Ph^P-Mx.<Csh'\\f\\x- (46) 

In other words, there holds 



PhA*-^ - A*-^Ph 



< Csh' . 

L{X,Xn) 



Remark 4.1 Compare to |LT06] . the important point to note here is the 
appearance of the function {A*)~'^^'2. in (35), (36) and (42). 

According to |LT06j . the inequality (22) make sense since 7 < | and thus 
imPh C D((-i*)i/2) c D{{-A*)^). 

In our context, on account of 7 > |, the inequality (36), which is sim- 
ilar to inequality (22) in |LT06j . only make sense if we add the functional 

(^j\*)-'y+2 in order that im{A* ^ Ph) C D{[—A*)^). The choice of the func- 
tion {A*)~^~^'2 seems to be the best adapted to our theory. 

Namely, we give here for instance about the functional [A*)~^^'^ through 
the heat equation with Dirichlet boundary control as follows 

y = /\y + c^y in (0, T) x Vt 

1/(0, .) =yoir).VL 

?/ = u in (0,T) X r = E. 

Set X = L'^{Q) and U = L^(T). It can be written in the form (4), where 
the self-adjoint operator A : D{A) — )■ X is defined by 

Ay = Ay + c^y : D{A) = [] H] ^ L\ 

In this case, the degree of unbounded of B is 7 = | + e(e > 0) (see jLTOOj . 
section 3.1). 

We may take A as follows 

Ah = -Ah, D{A) = H^f]Hl 
Therefore, (i*)-T+l = {A*)-^+' : H\n) Hl+'{n). 
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Remark 4.2 By means of the condition of the degree of unbounded of 
operator B and (33), we imply that B is admissible. 
Indeed, we have 



sup 



< B*S*{T - s)x,u{s) > ds 



< 



\B*sityijfdt)p 



^dt)pdt (P7 < 1) 



Remark 4.3. It is easily seen that assumptions (i^s) (except for the 
inequalities (35), (36), (42)) and hold for most of the classical nu- 

merical approximation schemes, such as Galerkin methods, centered finite 
difference schemes,... Additionally, by using some approximation properties 
and the properties of the functional {A*)~'^^^^'^, we prove that the inequal- 
ities (35), (36), (42) also hold for most of the above classical schemes (see 
the proof in Section 5). As noted in |LT00j . the assumption H41 of uniform 
analyticity is not standard, and has to be checked in each specific case. 

The main result of the article is the following : 

Theorem 3. Under the previous assumptions, if the control system y = Ay+ 
Bu is exactly null controllable in time T > 0, then there exist P > 0, hi > 0, 
and positive real numbers C, C satisfying 



IP 



< 



< 



'0 

C'\ 



\ldt+h^\\^,r^^ 



(47) 



for every h G (0, hi) and every iph ^ ^h, (1 < P < 2). 

In these conditions, for every yo G X , and every h G (0, hi), there exists 
a solution (fh G Xh minimizing the functional 

Jhi^h) = - [ \\Ble'^*^i;f,\\ldt+-h^\\i;hr^^+ < e^^^^^,P^i/o >x„ (l < P < 2) 
P Jo P 

(48) 

and the sequence {QhUh)o<h<hi , where the control Uh is defined by 
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for every t G [0,T] converges weakly (up to a subsequence), in the space 
L'^{0, T; U) to a control u such that the solution of : 

Ay + Bu, y{0) = yo, 



satisfies y{T) = 0. For every h e (0, hi), let yh{ ) denote the solution of 

Vh = ^hVh + BhUh, yh{0) = PhVo- 

Then, 

• yh{T) = -hf^ WifhW^'"^ (fh; 

• The sequence {Phyh)o<h<hi converges strongly (up to subsequence) in 
the space L'^{0,T; X), to y(.). 



Futhermore, there exists M > such that 

f \\u{t)\'^ < M^'^-'^ \\yo\i^-'\ 
Jo 

and, for every h e (0, hi), 

r\\unit)ru.<M^'^-'^ \\y4f'-'\ 
Jo 

\\yh{T)\\^, < M'l^^-'^h!"'^ \\yo\U^'^'^ . (49) 

Remark 4.4 The left hand side of (47) is uniform observabihty type 
inequahty for control system (2). This inequality is weaker than the uniform 
exact null controllability. No attempt has been made here to prove uniform 
exact null control for the approximation systems (2). 

Renicirk 4.5 A similar result holds if the control system (1) is exactly 
controllable in time T. However, due to assumption {Hi ), the semigroup 
S(t) enjoys in general regularity properties. Therefore, the solution y(.) of 
the control system may belong to a subspace of X, whatever the control u is. 
For instance, in the case of the heat equation with a Dirichlet or Neumann 
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boundary control, the solution is the smooth function of the state variable 
X, as soon as t > ,for every control and initial condition G L^. Hence, 
exact controllability does not hold in this case L^. 

Moreover, one may wonder under which assumptions the control u is the 
control, is defined by (27), such that y(T)=0. As in |LT06j . the following 
proposition give an answer: 

Proposition 1. With the notations of theorem, if the sequence of real num- 
bers llV^hllx^, < h < hi, is moreover hounded, then the control u is the 
unique control, is defined by (27), such that y(T)=0. Moreover, the sequence 
{QhUh)o<h<hi converges strongly (up to a sequence) in the space L'^{0,T;U) 
to the control u. 

A sufficient condition on yo G X , ensuring the boundedness of the se- 
quence {\\(ph\\xJo<h<hi, is the following : there exists rj > such that the 
control system y = Ay + Bu is exactly null controllable in time t, for every 
t E [T — ri,T + rj], and the trajectory t i— )■ S(t)yo in X, for t E [T — ri,T -\- rj], 
is not contained in a hyperplane of X. 

Other sufficient condition on control u, also ensuring the boundedness of 
the sequence (||v^/i||x;^)o</i</ii; is the following : there exists rj > such that 
the control system y = Ay-\-Bu is exactly null controllable in time t, for every 
t G [T — ?7,T + r/], and with the control u is defined as (21), the trajectory 
t\-^Sit — ^)Bu{C,) in X, for t E [T — ri,T+ri], every E (0,t) is not contained 
in a hyperplane of X. 

5 Proof of the main results 

1. The proof of theorem: 

Proof. For convenience, we first state the following useful approxima- 
tion lemma, whose proof readily follows that of |LT06] . [LTOO] . The 
proof of this lemma is provided in the Appendix. 

Lemma 3. There exists Cg > such that, for all t E (0, T] and h E 
(0, Hq), there holds 



(50) 




(51) 
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for every 6 e [0, 1] . 



(52) 



We carry out proving the theorem as follows: 

The degree of unboundedness 7 of the control operator B is lower than 
i, there exists 6 e (0, 1) such that < ^ + (1 - < i. 

For all h e (0, Hq) and -0^ G Xy^ we have 



dt 



u 



u 



+ 



B*S{tYPhiJ, 



dt. 



(53) 



We estimate two terms of right hand side of (53). 

The control system is exactly null controllable in time T, then there 
exists a positive real number C > such that 



dt> c 



X 



We have the following inequahty 

\yP-zP\ < p{yP-^ + zP-^) \y-z\, 



(54) 



(55) 



where y,z & R'^,p > 1- 

Indeed, we apply mean-value theorem for f{x) — x'^{p > l,x e Ef^) , 
there exists ^ e {y,z) such that 

ly'-z^l = \f{0\\y-z\ 

= p\C^-'^\-\y-z\ 



We apply the above inequality and make use of (40), (28), (44), (50) 
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to obtain 



< p{ PhS{TyPhiPn 



p-1 



Xh 



Xh ^ 

PhS{TyPniJh 
p-l^ 



Xh 



X,, 



< p{C,C^e-'\M\7h +C7Uhrx~ )- PhS{TrPf,,l;^-e^''^ij, 



Xh 



Therefore , from above estimate and (39), we get 

h^^'Mlxh-^^^^'W^i^^Wxh ^ PhS{TyPh^Ph 



Xh 



<ci 



s{TyPh^i 



X 

(56) 



Combine (54) with (56) we have: 



T 



B*S{tyPhiPh ' dt > Cislle^^^^^ll^^ - Cuh'WMxh- (57) 



For the first term on the right hand side of (53), one has, using (33), 
(51), (52) and applying the inequahty (55) 



< P{ 



B*S{tyPhiph 



u 



p-i 



u 



+ 



X 



B*s{tyPhiPh 



p-i 



u 



B*s{tyPhi^h 



u 



QhBle'^Hh - B*S{tyPnilJh 



C 



ip-l 



X 

^s(l-7)9 



U 
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We have 7 < ^(p > 1), therefore ^ + (1 — 9)^ + ^{p — 1) < 1 and we 
can get, by integration, 



Therefore, 

T 



u Jo 

(58 



We choose a real number /3 such that < /3 < s(l — 7)^. Combine 
results (53), (57), (58) we have inequality (47). 



For h e (0, /ii), the functional Jh is convex, and inequality (47), is 
coercive. Therefore, it admits a solution minimum at (ph G so that 

where GhiT) = \\Ble*^*hiph\\^''^e^'^''^BhB'^e^'^hdt is the Gramian of 
the semidiscrete system. 

With Uh{t) = ||5*e(^-*)^^(^ft||^~^S^e(^-*)^'^¥?;, is chosen then, the so- 
lution Uhi-) satisfies 



Jo 

Note that, since J^(0) = 0, there must hold, at the minimum, Jhi^Ph) < 
0. Hence, using the observability inequality (47) and the Cauchy-Schwarz 
inequality, one gets 

J 

and thus, 
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As a consequence, 



^0 



^ Ble^^l^4l^ < (^)V(--)(||P,yoirx?-^^), (60) 



and \\<Ph\\x^ < {^y^^~^W\Phyo\\x!!'~^^). and the estimates (49) fol- 
low. □ 

2. Proof of proposition 

Proof. If the sequence ( Phfh )o<h<hi is bounded then up to a sub- 

sequence, it converges weakly to an clement (p E X. It follows from the 
estimate (52) that u(t) = \\B*S{T - tyipf^"^ B*S{T - t)*ip for every 
t e [0,T]. Moreover, QhUh tends strongly to u in L^(0, T;t/). Indeed, 
for t G [0,T], 

QhUh{t) - u(t) 

- \\B*S{T - t)*ip\f~'^ B*S{T - t)*ip 
= ||S*eM^^(^,f {QuBle^^-'^^l - B*S{T - 1)*?^)^^ 
+ IIS^eM^^v'.f B*S{T - tYiPu^n - 

+B*S{T - t) V(||5*e(^-*)^^(/.^|f - \\B*S{T - t)*ipf-^). (61) 

Since the (ph are bounded, then the \\uh\\ are bounded. Prom that, we 
imply the ||S^e*^^~*)^ft(^^||^ are bounded . 

Using (52), the first term of right hand side of (61) tends to zero clearly. 
For the second term, for every t G [0,T] the operator B*S{T — t)* is 
compact, as a strongly limit of finite rank operators and since Ph<Ph — ^ 
tends to weakly to zero, it follows the second term of the right hand 
side (61) tends to zero. Furthermore, by applying inequality (55) we 
get 

||5*e('^-*)^'V/.|r"' - \\B*S{T - t)Vir-' 
< (2-p)(||i?*e(^-*)-^^v'jr"'+ ||5*5(T-t)V|r') 



x(||5*e(^-*)^^(^,-5*5(r-i)V 

As the ||5^e'^"^~*^^''93/i||'' are bounded and inequality (52) is used 
again. Hence, the third term tends to zero clearly. 
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The control u is such that y(T)=0, hence the vector must be solution 
of VJ*(^/') = 0, where J is defined as in Theorem 2. Since J is convex, 
(y9 is the minimum of J*, that is, u is the control such that y{T) = 0. 



We next prove, by contradiction, that the sufficient conditions pro- 
vided in the statement of the proposition implies that the sequence 
(llv'/i|lxh)o</i</ii is bounded. As the proof of the first sufficient condi- 
tion is found in |LT06j . we give here the proof only for the second 
sufficient condition. If the sequence {\\'fh\\xJo<h<hi is not bounded, 
then, up to subsequence, Ph{^h/ W^hWxJ converges weakly to $ in X, 
as h tends to 0. For every t e [T — r],T + rj], the control system is 
exactly null controllable in time t; and thus, from (60), the sequence 
Jg < B^e^'^'^^^'^ifh-, Qhu{s) >Uh ds is bounded, uniformly for h G (0, hi). 
Thus, passing to the limit, one gets 

t 

< $, S{t - s)Bu{s) >x ds = 0. 

This equality is equivalent to the fact that : there exists ^ G {0,t) such 
that < ^,S(t — ^)Bu{^) >x= 0. This contradicts the fact that the 
trajectory t i— )■ S{t — ^)Bu{C,), t E [T — ri,T + r]] and every ^ e (0, t), is 
not contained in a hyperplane of X. □ 

6 Numerical simulation for the heat equation 
with Dirichlet boundary control 

In this section, we give an example of a situation where the theorem 3.1 
are satisfied. 

Let Q C -R" be an open bounded domain with sufficiently smooth bound- 
ary r. We consider the Dirichlet mixed problem for the heat equation: 

y = Ay + c^y in (0, T) x Vt 

1/(0, .) = 1/0 in 
?/ = M in (0,T) X r = S, 

with boundary control u G L^(0,T; L^(r)) and y^ G LF'{VL). 

Set X = L^(fi) and U = L'^iT). We introduce the self-adjoint operator: 

Ah = Ah + c^h : D{A) = H'^n L^ip). 
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The adjoint B* e L{D{A*), U) of B is given by 

Moreover, the degree of unboundedness of B is 7 = | + e (e > 0) (see 
jLTOO] . section 3.1). 

1. One-dimensional Finite-Difference semi-discretized model: 

We next introduce a semi-discretized model of the above heat equation, 
using ID Finite-Difference. 

For simplicity, we set VL = (0, 1), F = {0, 1}, c=l and T=l. 

Given n G N we define h = > 0. We consider the following simplex 
mesh: 

= {xo = 0;Xi = ih, i = 1,...,N] Xn+i = 1} , 
which divides [0,1] into n+1 subintervals Ij = [xj,Xj^i] j = 0,...,n+l. Set 

Xh = {yeC%nH)}, 
U, = {yE C%T)}. 

Define i\ (resp., Qh) as the canonical embedding from into D({—AY^'^) 
(resp., from Uh to U). For Xh € Xh and Uh E Uh, set, Ph{xh) = 
Xh and Qh{uh) = uh- For y E Z)((-A)V2)' = H\ny, set P^y = 
{yi, ..,yi, ..,yn+2) where yi = y{{i - l)h) and, for u E U, set, QhU = 
..,Ui, .., Un+2) where Ui = u{{i - l)h). 

It is clear that the assumptions (-f^s.i) and (ifs.a) are here satisfied. Our 
aim is next to verify the inequalities in {H3 2) and (-^3.4). 

In order to get these inequalities, it will necessary to making use of the 
following usual approximation properties (see [LTOOj . section 5): 

{i)\\Uhy-y\\Hi(^n)<ch'~^\\y\\jjs^n),s<r + l, s - / > 0, < / < 1, 
and the inverse approximation properties 

(ii) \\yh\\Hc{n) < c/i-" \\yh\\l (fi), < a < 1. 

(iii) h-' \\y-Uhy\\^2^^ + \\{I -Uh)§;\\^,^^^ < ch'-^\y\\Hs^n)^ I <s< 
r + l,y E H'{n). 
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where r is the order of approximation (degree of polynomials) and 11^ 
is the orthogonal projection of L^(Q) onto Vh ■ 



First, by applying (i) we easily get the inequality (36) 

< Ch^ 



L2(0) 



i?2(Q) 



D(A*) 



< Ch'^WA*" 



in this case s=2. 



We next verify the inequality (37) as follows 



{{-Ay)(I-{A*) 'PnPnU 



< C 



{I - (A^f^hnPH)^ 



< Ch^ WH'WDiA* 

< Ch'^^-'^UlA* 



\X ' 



where we have used (i) with s=2, D{A*) = H'{n) and D{{-A*y) 



For the inequality (39), we apply (iii) with s=2 as 



(/ - QhQh)B*^ 



L2(r) 



(/ - QhQh 



dip 
dv 



L2(r) 



lD(A*) 

< Ch'^'-^^ \\A*^\\x . 
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For the inequality (43), by using (iv) and (40) we get 



B*{A* 



dv 



L2{r) 



< Ch-'^ 



< Ch—^ 



(i*) ^^''Phi^h 



L2(Q) 



Therefore, the inequahty (43) is satisfied for s=2, 7 = | + e. 



Moreover, the assumption {H4 2) is satisfied with s=2 (see |LTOO] ) . 
Hence, theorem 3.1 apphes, with /3 = 0.16, for instance. 



We then consider the following finite difference approximation of the 
above heat equation as follws 



Vj = J-^[yj+i + Vj-i - 2%] + c^Vj <t <T,j = 1, n 
%(0) = Vjo >J = 1> 



, n 



yo{t) = Vn+lit) = Uh ,0 <t <T, 

where y G , yo G -R"+^ , Uh e R and 



( ° 







0\ 





{c^h^ - 2) 1 










1 {c^h^ - 


2) 








{c^h^ - 


2) 


1 





1 














J 




/ 1 \ 







By 






VI/ 
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2. Numerical simulation 

The minimization procedure described in Theorem 3.1 has been imple- 
ment for d=l, by using a simple gradient method that has the advand- 
tages not to require the complex computations and this method can 
applied with any power p. However, the computation of gradient of Jh 
is very expensive since the gradient is related to the Gramian matrix. 



name 


Sh 


h 


Vo 


lD-10 


10 


10-1 


yo{x) = e-^' 


lD-100 


100 


10-2 


yo{x) = e-^' 


lD-500 


500 


2.10-=^ 


yo{x) = e-^' 



Table 1: Data for the one-dimensional heat equation. 



name 


yhWx 


h^(P = 0.16) 


\\yk{T)\\ 


lD-10 
ID- 100 
lD-500 


0.1690 
0.7960 
2.0570 


0.6814 
0.4779 
0.3699 


0.4775 
0.4565 
0.4273 



Table 2: Numerical results for one dimensional equation for ^ = 0.16. 



name 




/r'(J = 2) 




lD-10 

lD-100 
lD-500 


4.4266 

4.8933 
5.0956 


10-2 

10-^ 


0.0111 

1.3467e-004 
5.5178e-006 



Table 3: Numerical results for one dimensional equation for /3 = 2. 

Numerical simulation are carried out with a space-discretization step 
equal to 0.005, with the data of Table 1. The numerical results are 
provided in Table 2 for beta =0.16 and Table 3 for beta=2. 

The convergence of the method is very slow. Prom the result of The- 
orem 3, the final state VhiT) is equal to —h^ \\iph\\^~'^ iph in which iph 
is minimizer of Jh- We note that the maximum value for which the 
theorem asserts the convergence is very small. For such a small value 
of /3 (for instance (3 = 0.16), h'^ converges to very slowly. It follows 
that yh{T) converges very extremely slow. 

In practice, the unique minimizer (p^ of Jh is obtained through the sim- 
ple gradient method in which the step is equal to 0.01 and the error 
e = 10-2 is taken. With the small value /3 — 0.16, it took a long time 
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to achieve the results in Table 2. Namely, for n=500, it took more one 
month to get the result after 1000 iterations. It is clearly seen from Ta- 
ble 2 that the convergence of yh{T) is very slow. These results illustrate 
the difficult in using our method to compute controls. Although the 
case beta=2 is not covered by our main theorem, the method seems to 
converge for this value for beta and we provide hereafter the numerical 
results in Table 3. Since the value of the beta is greater, the convergence 
is quicker. 

Appendix: proof of lemma 

Proof. • First of all, we will prove (51) 
For every e D{A*), one has 



< 



u 



u 



+ 



B\A*) ''^'^PhPhSitYiP 



u 
(62) 



We estimate each term of the right hand side of (62). Since Bl 
QhB*{A*) ^Ph and thus, using (40) (42) (44) one gets 



QnBle'^np^^ 



u 



On the other hand, from (28), (40), (42) one obtains 



(63) 



B*{A*)~'^'^PnPhS{tyij 



u 



< Ch-^^WP^Sitmi^^ 

< C5Ceh-^'\\S{ty^\\^ 



(64) 



Hence, combine (63), (64) with (62), there exists Cio > such that 



QnB*e'^tp^^-B*{A*) "'P^P^SitYi^ 



7+ 



(65) 



u 
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for every ij; e D{A*), every t e [0, T], and every h e (0, ho). 

Moreover, we get another estimate of this term. By using (33), (38), 
(40) , (42), (50) one gets 



^ tAi ^ . ^.^-7+^ ^ 



u 



< 



'1+ 



QhQhB*{A*) ■'Pn{e'''^Pu^-PnS{ty^) 



u 



+ 
+ 



QuQHB*{{A*r^^' PhPh - I)S{trij 



{QnQh - I)B*S{ty^ 



u 



B*{I-{A*)~^^~'PkPh)S{ty^ 



< C,Ceh'''\\e'^*'^Phip-PhSity 



u 



+C5C3 



-Ay {{A*) 'PuPh-i)S{ty^i^ 



+W(i-^)||A*5(i)*V|lx 



+C3 



-Ay {{A*) 'PhPh-i)s{tyi/j 

/i«(l-7) 



< c 



11" 



1 



X + (^3(C^5 + 1) + i)c,h<^-''^A*s{tynx 

(66) 



Then, raising (65) to the power 1— 7, (66) to power to 7 and multiplying 
both result estimates, we obtain 



QhBle'^'^Ph^ - B*{A*)~'^~'P„PhS{ty^ 
Hence, 

Qf.Ble'^-'^Pn^ <%||^||^ + 



< 



C 



u 



12 



B*{A*) '^~'PnPhS{tyi: 



■ (67) 



u 
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From (29), (33), (36) one yieds 



< 



u 



+ ||5*^(t)>|| 



u 



< c 



13- 



(68) 



Combine (67) with (68) and by setting ip = Phiph we get (51). 

Finally, we prove (52). On the one hand, reasoning as above for obtain- 
ing (66), we get 



/,s(l-7) 

< c- 

u - t 



XI 



(69) 



for every ip G D{A*), every t G [0, T] and every h G (0, ho). 
On the other hand, from (51) and setting ip = Phiph one obtains 



QhBle'^*^Ph^lJ-B*S{ty^lj 



< 



QhBle'^*^Pn^lj 



< ^ 11^.11 +C3 

c .. .. 

< — 



+ \\B*s{ty 
-A*ys{tyij 



^hWx 



X 

(70) 



Raising (69) to the power 6*, (70) to the power 1 — and multiplying 
both resulting estimates, we obtain (52). 

The proof of the inequality (50) is found in |LT06j . [LTOOj . 

□ 



7 Conclusion 

We have shown that the appropriate duality techniques can be applied to 
solve (3), namely the Fenchel-Rockafellar theorem. 

Additionally, it is also stated that under standard assumptions on the 
discretization process, for an exactly null controllable linear control system, 
if the semigroup of approximating system is uniformly analytic, and if the 
degree of unboundedness of the control operator is greater than \ then the 
unform observability type inequality is proved. Consequently, a minimiza- 
tion procedure was provided to build the aproximation controls. This is im- 
plemented in the case of the one dimensional heat equation with Dirichlet 
boundary control. 
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Note that, we only stress our problem on the case 7 > 1/2. Some relevant 
problems for which 7 < 1/2 that are refered to |LT06j . 

One open question is given: how the above results change if we remove 
the assumption of uniform analyticity of the discretized semigroup. 
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